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Abstract 

The color-flavor transformation, an identity that connects two integrals, each of 
which is over one of a dual pair of Lie groups acting in the fermionic Fock space, is 
extended to the case of the special unitary group. Using this extension, a toy model 
of lattice QCD is studied: Nj species of spinless fermions interacting with strongly 
coupled SU(Ac) lattice gauge fields in 1 + 1 dimensions. The color-flavor transformed 
theory is expressed in terms of gauge singlets, the meson fields, organized into sectors 
distinguished by the distribution of baryonic flux. A comprehensive analytical and 
numerical search is made for saddle-point configurations of the meson fields, with 
various topological charges, in the vacuum and single-baryon sectors. Two definitions 
of the static baryon on the square lattice, straight and zigzag, are investigated. The 
masses of the baryonic states are estimated using the saddle-point approximation for 
large Nc- 
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1 Introduction 


In quantum chromodynamics (QCD), a hierarchy of scales is provided by ~ 1 GeV, the 
scale of chiral symmetry breaking, and Aqcd ~ 0.18 GeV, defined as the location of the 
Landau pole of the one-loop beta function. The running coupling constant increases from 
weak to strong coupling as the momentum scale is lowered from the perturbative regime 
above A^ down to Aqcd- 

In the past two decades a great deal was learned about the non-perturbative structure 
of QGD at scales between A^ and Aqcd- The guiding idea was to construct low-energy 
effective theories which encode the symmetries of the fundamental QGD Lagrangian. To 
obtain these effective theories, one may start from full QGD, and integrate out the high- 
energy degrees of freedom (quarks and gluons) in order to produce a low-energy effective 
action in terms of mesons and baryons. In this way it was possible to recover the chiral 
Lagrangian [Q, |, |] that had been introduced phenomenologically by Weinberg |^. 

In a more recent development, it was shown how to extract the effective long¬ 
distance degrees of freedom by starting from the lattice [0, §] formulation of QGD. In that 
approach it is assumed that the long-distance physics of lattice QGD (LQGD) can be 
described by a strongly coupled lattice theory. From the latter, one gets the continuum 
chiral Lagrangian by expanding the effective action in powers of the lattice spacing and 
external momenta. All the terms of the Gasser-Leutwyler continuum effective Lagrangian 
1^ can be recovered in this way |TO|. The lattice formulation is, however, deficient in one 
respect: by the technical difficulties with chiral symmetry for lattice fermions, the chiral 
anomaly is lost, i.e. for Nj massless quark flavors the chiral symmetry of the lattice effective 
theory is U(V/) rather than SU(V/). 

This type of approach was initiated in ||ll|, |l^ ; it relied on a “bosonization” of the 
strong-coupling LQGD action, and a large-Vc or large-dimension expansion. Technically, 
the heart of the method is the computation of integrals over the group SU(Vc) with Haar 
measure, weighted by Some general results for such integrals have recently been 

reviewed in |]T^ . 

A few years ago, an alternative kind of bosonization scheme was introduced [|^], relying 
on a mathematical formalism later called the “color-flavor transformation” [BAII. This 


transformation relates two different formulations of a certain class of theories. In condensed 
matter theory, the transformation has found a number of applications, among others to 
the random flux model 
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The color-flavor transformation in its original version applies to the gauge group U(Vc). 
For this group, all gauge singlets are of “mesonic” (or quark-antiquark) type. In order for 
baryons to appear, one needs to replace U(Vc) by the special unitary group SU(Vc). In 
Section ^ of the present paper we extend the color-flavor transformation to SU(Vc), by 
decomposing the (colorless) flavor sector of Fock space into disconnected subsectors labeled 
by the baryonic charge. 

In Sections tB we apply the formalism to a toy model of LQGD: Nf species of spinless 
fermions interacting with strongly coupled SU(W) lattice gauge fields in 1 + 1 dimensions. 
The color-flavor transformation yields a dual representation of this non-Abelian model. 
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Combining numerical computations with analytical considerations, we conduct a compre¬ 
hensive search for saddle-point conhgurations in various baryonic sectors with different 
topological properties. We use these conhgurations (without huctuation corrections) to 
estimate the mass of a single baryon in our model. In doing so we ignore the Mermin- 
Wagner-Coleman theorem (asserting that spontaneous breaking of continuous global sym¬ 
metries does not occur in 1 -I- 1 dimensions), by assuming the pattern of chiral symmetry 
breaking that is known to occur in the physical case of 3 -|- 1 dimensions. 

After the present work had been completed, we learned that the SU(iV) generalization 
of the color-havor transformation has also been worked out by Schlittgen and Wettig m- 


2 Color—Flavor Transformation for S\]{Nc) 


2.1 Group action on fermionic Fock space 

In this section we set up some algebraic structures, which are needed to establish the 
“color-havor” transformation for the special unitary group. Our discussion follows the line 
of reasoning of Ref. [|^] but is somewhat simpler, as we do not need the superalgebraic 
framework employed there. 

We start by considering a set of fermionic creation and annihilation operators and 
f\, which obey the canonical anticommutation relations 


{/t /b = 0. {/i. /b = 0. {/i. /b = sabs » . 


The lower index takes the values +a or —a, with range a = 1,... ,Nf, and the upper 
index takes the values i = 1,, N^. Having QCD in mind, we interpret the operators 
and as creation operators for “quarks” and “antiquarks” respectively; the index i 
corresponds to the gauge (or color) degrees of freedom and the index a labels the diherent 
quark havors. (The quarks are regarded here as being spinless.) The operators and f\ 
act on a Fock space with vacuum |0) and its conjugate (0|, by /)^|0) = 0 and {0\f^ = 0 for 
all A and i. 

We next consider the set of quadratic operators dehned by 

pij ^ ri rj rpij ^ ri p 

'‘~^+a,+b J+aJ+b ! -‘-^+a,—b J+aJ—b i 

_ p p Tpp _ fi p 

-^-a,+b ~ J-aJ+b ! ^-a-b ~ J-aJ-b ' 

The C-linear span of these operators has the structure of a complex Lie algebra, 0. More 
precisely, the operators obey the commutation relations of a set of canonical generators 

of the Lie algebra gl{2NfNc): 

lA'fl. JSSdI = Si%cE% - S“SoBEi’B . 

Thus we have a Lie algebra isomorphism from g[{2NfNc) (i.e. the space of complex matrices 
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of size 2NfNc x 2NfNc, with the Lie bracket given by the commutator) to 0: 


t : 0[(2iV^iV,) 
m 




tm ■= ^AbEaB ■ 

ij,AB 


This isomorphism lifts to an isomorphism of the corresponding complex groups: 

T : GL{2NfN^) G , 


M = exp(m) 


Tm = exp(t^) , 


( 1 ) 


which forms a (reducible) representation of G\j{2N on Fock space. The representation 
is single-valued (which means there are no U(l) obstructions from the multi-valuedness of 
the logarithm) as the spectrum of each operator is the set {0,1}. 

The Lie algebra Qi{2NfNc) has two subalgebras 5 l(A^);) and Q[{2Nf) which are embedded 
in a natural way: a matrix X G qI{Nc) is identihed with ^Nf G)X, and a matrix Y G g[(2iVj) 
with Y ® ^Nc- Through these embeddings, qI{Nc) and Ql{2Nf) form a pair of maximal 
commuting subalgebras of 0[(2iVjiVc), also known as a “dual pair” |^. The subgroups 
GL(A'"c) and GL(2A^j) are embedded into GL(2A^jiVc) in the same way. Their adjoint action 
on the fermionic creation and annihilation operators is described in Appendix 

We dehne the color group to be the subgroup SU(Ac) of GL(Ac), and the flavor group 
to be the subgroup Y{2Nf) of GL{2Nf). GL(Wc) contains an extra U(l) subgroup which 
lies outside the color group and, being generated by the unit matrix, commutes with the 
whole group GL{2NfNc). This U(l) is generated hy Q + Nf where 




AA- Nf - Jj- y^Xf+afla 

^ a.i 


f-af-a) 


( 2 ) 


counts the difference between the number of particles and antiparticles: Q = ^(W+ —W_). 
In contrast, the operator giving the total number of particles, N = 'YlaXf+af+a + f-af-a) y 
does not commute with the generators of gl{Nf). We will call Q the baryon charge operator. 


2.2 Prom color group integrals to flavor group integrals 

Let and be two independent sets of Grassmann variables, referred to as “quark 
helds”, and consider the color group integral 

JsV{Nc) 

The Haar measure dU of SU(Ac) is understood to be normalized by /gu(jv) ~ 

also adopt the convention that repeated occurrence of an index implies summation. 

The color-flavor transformation will replace the integral (^ by an integral over the 
flavor group Y{2Nf) . A key step in doing the transformation is to interpret as 
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the matrix element of an operator V that projects on the colorless sector (or flavor sector) 
of Fock space. This sector is the subspace of all states |flavor) which are invariant under 
the color group: Tt/|flavor) = |flavor) for all U G SU(iVc). 

The hrst step towards the color-flavor transformation is to express the projector V as 

V = [ dUTu. (4) 

JsV{Nc) 

Let us now introduce the fermion coherent states 

(Tl = (0| + V'+JU , |T) = exp(/l„V’-a + MlaM ■ (5) 

By making use of the hrst set of relations in Appendix A, it is straightforward to show 
that 


(f |Tt;|T) = exp(^;,f/*^'^i„ + 

for U G SU(A'"c). This yields the simple formula 

Z(V^,^) = (T|P|T). (6) 

To express as an integral over the havor group, we will derive an alternative 

representation of the projector V, as an integral over coherent states of the havor sector. 


2.3 The flavor sector 

The subspace of states in Fock space which are invariant under U{Nc) was described in 


r^. It consists of the vacuum and of mesonic excitations on top of it. The prototype of 


such an excitation is the “one-meson” state 


|m.») = Ka.-M = E /; jy |0) 


By the multiple action of the gi{2Nf) generators (where we have gone back to using 

the summation convention), one can build states containing up to N^Nf mesons, with 
diherent havors. These states are automatically U(A'),)-invariant; conversely, all U(A'"c)~ 
invariant states are linear combinations of such multi-meson states. The group \J{2Nf) 
acts irreducibly on this invariant subspace. 

The set of SU(A'),)-invariant states is larger. To obtain it, one relaxes the constraint 
Q\‘ip) = 0. Thus there exist colorless sectors of Fock space on which the central generator Q 
takes a non-zero value. These sectors contain the baryons, which are totally antisymmetric 
combinations of Nc quarks. A baryon with flavors ai,... , is dehned as 




_;r. . fU 

jY j^n...*iVcJAi 


•ejo), 


(7) 
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where the Ak = are taken either all positive (baryon), or all negative (antibaryon). A 
matrix g G GL(A'"c) acts on this state simply by multiplication with Det( 5 ') (resp. 

Therefore, the state is invariant under the color group SU(Ac). 

The above baryon (resp. antibaryon) is an eigenstate of the baryon charge operator Q 
with eigenvalue +1 (resp. —1). Acting on it with the generators E'^2_b flavor algebra 

0 [( 2 A^j), one builds other colorless states with the same baryon number, which form an 
irreducible subspace for U( 2 A^j-): the one-baryon (resp. one-antibaryon) sector. 

The one-baryon sector can be generated from the state (i with all ttj = 1. One can 
similarly build Q-baryon (resp. Q-antibaryon) states from 


Q 


Q 


|5o) = | 0 ), |i?_Q) = n/- 


/-alO) 


( 8 ) 


a=l 


a=l 


The values of the baryon charge range from —Nf to Nf, according to Pauli’s exclusion 
principle. As with Q = ± 1 , acting on \Bq) with the algebra g[( 2 A^j) builds the full Q- 
baryon part of the flavor sector, so the group U{2Nf) acts irreducibly on this part. This 
can be proved by using the dual-pair property of the subalgebras 0 [( 2 A'j-) and qI{Nc), as 
exposed in |T^ . 


To summarize, the flavor sector of Fock space decomposes into 2Nf + 1 subsectors, 
characterized by their baryon charges Q. Each sector carries an irreducible unitary repre¬ 
sentation of the flavor group U( 2 A^/). 


2.4 Coherent states 


Having decomposed the flavor sector as described above, we can now express the projector 
P in a different way. For this purpose we will use coherent states, in the spirit of Perelomov 
IP| . On each subsector with a hxed baryon charge Q, we consider the generalized coherent 
states built by the action of G = U(2A^/) on the reference state \Bq), i.e. the states 


Vj € G, VQ = -N,, | 9 g) S T,\Bq), (jg |S (Bg|Tt . (9) 

The crucial property of coherent states we will now use, is that they supply a resolution of 
unity. Because of the irreducibility of the U(2A^/) action on each Q-subsector, the operator 


'PQ = aQ [ dg\gQ){gQ\ ( 10 ) 

Jg 

coincides with the orthogonal projector on that subsector, the only provision being that 
the normalization constant aq be chosen appropriately. Indeed, the operator Vq trivially 
commutes with every element of the flavor group; Schur’s lemma then ensures that it 
is proportional to the identity on each irreducible space of this group, therefore on each 
subsector with hxed baryonic charge. Owing to orthogonality, Vq vanishes on all subsectors 
with Q' 7 ^ Q, whereas it is the identity on the Q-subsector if we take 
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Some particular values of the constant (namely oq, a±i) are computed in Appendix 
For the matrix element (^) of the projector V on the full flavor sector, 

Nf 

© -Pq. 

Q=-Nf 

we now have a new representation: 


Nf 


X] ttQ / dg{'^\gQ){gQ\^) . 

r\— AT. G 


( 12 ) 


0—IV, 


To compute the overlaps {'^Igg) and (^fglT), it is convenient to use a Gauss decomposition 
of G = U(2A^/): almost any matrix g = G G can be factored as 


A B 
C D 


1 Z\ (A 0\ (I 0 

0 1j I 0 D [z 1 


(13) 


where the relations Z = BD~^, Z = D~^C, and A = A — BD~^C hold. The decomposition 
becomes singular if D does, but this happens only on a submanifold of codimension one 
(and hence measure zero) of G. The unitarity of g implies Z = —D'^Z'^A^ and allows to 
write the central matrix in the form 


[A 0\ _ /(l + ZZt)i/2 0 \ /W 0\ 

VO 0 (i + ztz)-i/7 VO Vj 


(14) 


U and V are unitary, so is an element of the diagonal U(A^/) x U(A^/) subgroup of 

G, which we call H. It can thus be shown that the elements g of an open dense subset of 
G are in one-to-one correspondence with the triplets (Z, G, V), where the pair diag(G, V) 
is an element of FT, while Z represents a point in the coset space G/H and can be any 
complex Nf X Nf matrix. Moreover, the Haar measure dg of G factorizes as 

[ dg= [ d{gH) [ dh = [ dgi{Z,ZA [ dUdV. (15) 

JG Jg/H Jh Jh 

Both did and dV are normalized Haar measures on U(A^/), and 
d/i(Z, ZA = Gnj Det(l + ZZA~^^f JJ 


is the normalized invariant measure on G/H. The normalization factor G^f is computed 
in Appendix see Eq. (^) . 
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We now explain how to use this decomposition to compute the overlaps. The Gauss 
decomposition ([I3D carries over to any representation of G, so for every g & G we can write 
the operator Tg as 


^diag(d,D) 


(16) 


where ( 


1 Z 
0 1 


and ( 


1 0 
Z 1 


According to the relations given in Appendix 


the factors and act trivially on the reference states: 


yQ = -N,,... ,N,: T{|Bo) = |Bo>, = {Bq\ . (17) 

The action of the block-diagonal operator is slightly more subtle. Using the third set of 
relations given in Appendix A, we get 

rd,.,(i,D)| 0 > = (Det£>)"-| 0 ) , 

^diag{A,n))l-®l) = (DetD)«"n'4«i/hlO>. 

i=l 

No 

= (DetD)""n(-D'‘)i«/iJO) ■ 


(To make sense of these formulas one must remember that we are using the summation 
convention: the flavor index a under the product is understood to be summed over.) These 
formulas directly yield the desired overlaps with (T|: 


Nc 

(T|^o) = {Bet , 

i=l 

No 

(T|^i) = {BetD)^'^Yl'^+cNiew{'4’XaZab'il^-b) i 

i=l 

No 

i=l 

as well as the overlaps with |T): 

Nc 

(^ol^) = {Bet abi^Xb) 1 

i=l 

No 

(^i|T) = {Bet D^)^^Yl^ic'^+c^^v{'^-aZKb'ip\b) , 

i=l 

No 

(9_.|>f) = (DetDt)'''n''’-o(r>-‘)Uxp(V-yzt„t^y) . 
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The overlaps with the coherent states Iqq) containing more than one baryon (|Q| > 1) can 
be computed in the same way; in front of the exponential factors, there will be IQI similar 
products, with flavor indices 1 ,... , IQI- 

We now insert the above expressions for the overlaps into ( 0 . and use the factorization 
(p!5|) to arrive at an integral over triples {Z,hl, V). Leaving the Z-integral for later, we next 
carry out the integrations over the unitary matrices U and V. They enter in the overlaps 
via the matrix elements of A and D] see Eq. ( 0 ). To simplify the notation, we hrst perform 
a flavor rotation on the Grassmann helds: 

= ^U^i + zz%, , = (Vi + z^z),J\. 

The integrals we need to compute then read as follows (assuming Q > 0): 


Q Nc 


Xq ((/)+, (/>+) = ag dUYl W{(t)\aUac) {(l>\bKb) . 


U(iV/) 


c=l i=l 


X-q( 0 -, 0 -) = «- 


Q 


V{Nf) 


Q N, 

dvlllli^\v;:){ci>UVbc). 

c=l i=l 


(18) 

(19) 


We also set Xo = ao, and Xq(V’, i’\ Z) = Xq(<I>+, <t>+), and X-qO’, Z) = X-Q{<t>-, <!>-)■ The 
function xi{.4>+, </>+) will play a distinguished role in the lattice gauge theory application in 
Section and we therefore evaluate it explicitly in the next subsection. 

The integrations over H having been done, we are left with an integral over G/H, 
i.e. over a Z-dependent integrand, in each Q-subsector. Putting everything together, we 
hnally arrive at the following identity: 


JSV{N,) 

Nf 

E 

Q=-Nf 




’c f 


dn{Z, Z^) xq(V',' 0 ;Z) 


exp{^lJ\^Zab^/Jlb + i^lbZ\ai^ia) 


Det(l + ZZt) 


N, 


( 20 ) 


which is called the color-flavor transformation for SU(iVc )5 and is the central result of the 
present section. Note that the right-hand side of the transformation has the attractive 
feature of organizing the contributions according to the different baryonic sectors. 

An effective action in the bosonic variable Z can be obtained by doing the (Gaussian) 
integral over the Grassmann fields. This will be done in a lattice gauge context in Section 







2.5 Evaluation of xi 

In this subsection we evaluate the coefficient 

Xi(0+,0+) = tti / dU 

Jv{Nf) 

Only the hrst column of the unitary matrix U occurs in the integrand, so the integral 
is effectively over a unit sphere in iV/-dimensional complex space, 

Parametrizing the latter by a complex vector z = {zi,... , z^f) with unit norm \z\ = 1, we 
have 




Xl( 0 +, 0 +) = « 1 - 


n 

2 = 1 


dn{z, z) n (</>+a^a)( 0 Vfe) 




where dfl{z,z) is a U(A^j)-invariant measure on the unit sphere \z\ = 1. By homogeneity 
in ^ and z, we may use the trick of replacing the numerator and denominator by integrals 
over , with a Gaussian weight function e“l^l included in the integrands. The answer 
then easily follows from Wick’s theorem: 


Xi{(f>+,(f>+) = ai 


= «! 


(Nf-iy. 
{N, + Nf-1)\ 

(Nf - 1 )! 

(iV, + iV;-l)! 


Nc 


sgnaj] 






i=l 


^ sgna 


(b 

+6 ’ 


( 21 ) 




2 = 1 


where ©at^ denotes the group of permutations of the numbers 1 ,... , W 


3 Color—flavor transformation on the lattice 


We consider a Euclidean SU(iVc) gauge theory in 1 + d dimensions placed on a hypercubic 
lattice with lattice constant a. The fermions 'ipl{n), with colors i = 1,... , W and flavors 
b = 1,... ,Nf, are put on lattice sites labeled by n = (no,... , n^), while the gauge matrix 
variables t/(n+|) = exp (^iagA^{na + ^)) G SU(W) are placed on lattice links n + /i/2 
(we label links by their middle points)^ starting from sites n in any of the directions 
/i = 0,... , d. In the limit of a strong gauge coupling //, the gauge theory has the partition 
sum Z = j Y\di>{n)dyj{n)Z{tlj with 




'SU(7V^) 


dU{n+^) 


0*^ C/,'0,’0 


(n+f) 


( 22 ) 
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The fermions on two neighboring sites n and n + /t are coupled through the gauge helds 
on the connecting link n + /i/2 in a gauge-invariant way: 

= Y + h) - + /i)17^-^*(n+f)i/fe(n)) , (23) 

while the (bare) quark mass m couples the fermions diagonally 

5'm,v.,^(’^) = a‘^’^^mi/(n)'0(u) . (24) 

We are not going to worry here about the fermion doubling problem and will restrict our 
considerations to this naive discretization of the fermionic action. Also, for simplicity we 
do not take into account the spin degrees of freedom, leaving their inclusion for a future 
publication. 

We rescale the fermionic helds so as to absorb the prefactor 0*^/2. This just adds a 
global prefactor to Z, and has no effect on the physical quantities. The SU(W)Aiitegral 
over U on each link is then identical to (j^) after the following substitutions: 

?/+='0(u), ='ip{n + fi), '4>_ = 'rp{n), -0-='^(’^ + A)- 


On each link, we perform the color-havor transformation (pO]) , thereby introducing a 
complex “havor matrix held” Z(n+|), Z\n+^). The outcome of the transformation reads 


{<9} « 


2am'i/>(n)'i/;(n) 


n 




dll(Z,Z^(n+l))x‘i 




Z,'iIj, i}j V 

) Det(l + ZtZ(n+f))«- 


(25) 


where the sum on the right-hand side extends over all possible distributions {Q} of baryonic 
charge (actually, baryonic hux) over the links of the lattice. The color-havor transformed 
action on a link n + fi/2 is 

‘^z,y,y(«+t) = i’l{n)Zab{n+^)ijl{n) + V>^(u + + fi) , (26) 

and the y-coefhcients are = cto and (with Q = Q{n+^) > 0) 

= XQ(n+|) + ZZ'f{n+!^), ^Jl + ZZ^{n+!^)ll){n + /i)^ , (27) 

= X_Q(n+f) (^\/1 + Z^Z{n+^)'il;{n),i!{n + il)^l + ZtZ(n+f)^ . (28) 


The fermions are now coupled through their flavor indices, whereas in the original action 
the coupling had been mediated by the color degrees of freedom. Moreover, the coupling 
has become ultralocal: the fermions at a site n couple only to one another, via Z{n+^), and 
so do the fermions at site n + fi, via Zl(n+|). Correlations between neighbors are solely 


10 







due to the relation between Z and by Hermitian conjugation. A graphical description 
of the change of coupling scheme is given in Fig. |^. 

The partition function (^) is a sum over all conhgurations of baryonic fluxes )}• 

For most of these conhgurations, the Grassmann integral vanishes identically. To see that, 
we expand the integrand for a given conhguration into a polynomial in the Grassmann 
helds, and count (for each site n) the number of fermions '0(n), '^(n) in the various mono¬ 
mials: 


• For every direction /x, the coefficient x‘^(n+§) contains Nc\Q\ Grassmann variables 
"0* (n) if Q > 0, and the same number of Grassmann variables 'ipa{n) if Q < 0. 

• For the coefficients fhe situation is the same, except that 'ip{n) and 

switch roles. 

• Each term of the expansion of involves as many '0(n) as 

The Grassmann integral f d'ip{n)d'ip{n) extracts the coefficient of the top-monomial. 


/ 


d%lj{n)d%lj{n) 


N, Nj 

nif''^ 

2=1 a=l 


[nWa{n) 


= 1 


setting all others to zero. This monomial contains as many 'ip{n) as 'ipin). Hence, in view 
of the counting above, the contribution from a conhguration {Q{n+^)} vanishes unless the 
following condition is met: 


d d 

'^Q{n+^) = '^Q{n-^) . (29) 

fi=0 fi=i} 

The physical meaning of this equation is conservation of the baryon current: the (algebraic) 
number of baryons “arriving” at the site n (from the links n —/i/2) must equal the number 
of baryons “leaving” the site (via the links n + fi/2). 

The general structure of the partition function (^) corresponds to the hadronic corre¬ 
lation function written in terms of colorless A"c-quark currents p0|, 

3.1 Integration over the fermions 

Based on the general considerations above, we perform the integration over the fermions 
sector by sector, and present below two particular cases: the “vacuum”, i.e. the sector 
where the baryonic hux Q(n+f) vanishes on every link n + fi/2, and a toy model of a static 
baryon on a mesonic background (with Q(n+|) = 1 along a time axis). 

In each case, integration over the Grassmann variables yields a purely bosonic effective 
action, which depends on the conhguration of the helds Z and ZF After computing 
these ehective actions, we will look for their saddle-point conhgurations to estimate the 
respective partition functions. 
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Before computing the effective actions in particular cases, we emphasize the conse¬ 
quences of chiral symmetry, which emerges in the limit of zero quark mass. The hyper- 
cubic lattice is bipartite, so it can be split into two sublattices according to the parity of 
1^1 Given this splitting, the effective actions S{Z, Zl) of all sectors are invariant 

under the following global transformation: 

\n\ even : Z(n+|) i— ^ UiZ{n+^)U 2 , Z'*'(n+|) i— UlZ'^{n+^)Ul , (30) 

\n\ odd : Z(n+f) ^ UlZ{n+&)Ul , z\n+^) ^ UiZ\n+&)U 2 , 

for any pair {Ui, U 2 ) G U(A^/) x U(A^/). Therefore, in each sector, the saddle-point conhg- 
urations in the chiral limit form a continous set (namely an “orbit”) generated by acting 
with the chiral symmetry group U(iV/) x U(A^/). As soon as the quark masses are turned 
on, this degeneracy disappears, and the saddle points become isolated. Equations (^0]) 
show that the helds Z, Zl(n+|) transform differently according to the parity of \n\. To 
stress this difference, we give different names to the helds on different sublattices: the helds 

living on the “even” lattice links will be called V, 1/1 (n+f), while the helds on the “odd” 

links will be denoted by IT, lTl(n+f). 

3.1.1 Vacuum action 

For the vacuum sector we have zero baryonic hux (Q = 0) everywhere on the lattice; the 
integral over the fermions, being Gaussian, is then easily done and yields 

J d'ip{n)d'4){n) -0) 

f d^i{Z,Z\n+^))} exp(-Vc5'vacuum[^]) , (31) 

where the result of the integration has been sent back to the exponent. The integration 
measure in curly brackets will be denoted by 'Z1(Z, Z^) in the following. The factor A), 
in the exponent comes from the color content of the fermions: since the action S'^('0,'0) 
does not couple fermions with diherent colors, the Grassmann integral is a product of A), 
identical integrals. The ehective action is 

d 

TrlnV(n+|) , (32) 

n n /i=0 

where 

d 

M{n) 2am + ^ (Z(n+§)Zl(n-|)) , 

1 - 1=0 

iV(„+f) = l + Z(„+|)Zt(„+|), 


(33) 

(34) 


Z 


vacuum 
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3.1.2 Static baryon action 

By the static baryon we mean the following distribution of baryonic fluxes over the lattice: 

1 along the links of the “world line” (or “string”) n = (f, 0,.. . , 0) G 
/i = 0, with f = 0,... , T — 1; on all other links Q = 0. This distribution satishes the 
current conservation law (p9D at all sites but the ends f = 0 and f = T of the world line. 
There it does, too, if we impose periodic (or antiperiodic) boundary conditions on the 
Grassmann helds: for a lattice of size T in the time direction, we set 'ip{T + 1, r) = "0(1, r), 
'ip{T + l,f^ = "0(1, r). We again write the partition function in the form (^I]), 

-2^baryon j T>?/>T>^/^Zbaryon(V’, V’) = j exp(-iVc5'baryon[2']) . (35) 

The effective action S'baryon contains the “sea” term Svacuuml-Z^], plus an extra part coming 
from the factors Xi along the world line of the baryon. These factors depend on the values 
of the Z field along this line and on the adjacent links, through the following matrix: 

G = Ar(io)M(io)-^iV(f6) ■ ■ ■ A^((r-|)6)M((T-i)6)-^iV((T-|)o)M(T6)-^ . (36) 


(We use the abbreviation fO = 0 + fO to denote the sites or links on the world line of 
the baryon.) This product of matrices runs over all sites n on the baryon world line (it 
is expressed as a “quark propagator” along that line). In Appendix y we show that the 
effective action takes the form 


A^C Sbaryon [Z] _ j ^ f + Nf 1 




Nr, 


-T 


Nr 


Y, m) n (Ti'G') 

^eSiVc 


l\Cl{a) ^-NrS. 


c ^vacuum 


[Z] 


1=1 


(37) 


In the non-vacuum factor of (|37|) , a runs over all conjugacy classes of the group of 
permutations of the set {!,..., W}- Every representative of the class a can be decomposed 
as a product of cycles of various lengths /, such that q((t) cycles of length I occur; thus, 
each class a is uniquely specified by the sequence {q}, or equivalently by a Young diagram. 
The weight factor A/'((t) is simply the cardinality of the class d, and is given by 

N ' 

A7(d) = -. (38) 

n 

1=1 


For the lowest numbers of colors the explicit expressions are 


Nr 1 . 5'baryon ‘S'vacuum lu TrG T COUSt , 

Nr = 2 : S'baryon = S'vacuum “ ^ ((TlG)^ + TrG^) + COUSt , (39) 

Nr = 3 : Sbaryon = S'vacuum ~ ^ ((TrG)^ + 3TrG^TrG + 2TrG^) + COUSt . 
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The constants, which are not given above, make contribntions to the baryon mass, so they 
need to be taken into acconnt in the hnal answer. 

In the following section, we look for the saddle-point conhgnrations of the effective 
actions 5'vacuum cllld 5'baryon- 

3.2 Saddle—point equations 

In the two sectors that we are interested in - the vacnnm and the static baryon - we wish to 
compnte, or at least estimate, the partition fnnctions (^) resp. (|35D . Since we are nnable 
to provide an exact answer, we will treat both integrals in a saddle-point approximation, 
valid in the limit of a large nnmber of colors Nc- For both the vacnnm and the static 
baryon, we will restrict onrselves to a pnrely classical approximation, which is to say we 
will identify the saddle points, evalnate the action fnnctional on them, and approximate 
the partition fnnction as ^ ~ 0 -S(-^s.p.) xhns we neglect all loop corrections, which are of 
higher order in 1 /iVc. 

In the vacnnm sector, where Nc appears explicitly as a factor of Svacuuml^], the saddle- 
point approximation is fnlly jnstihed in the large-iVc limit. The sitnation is less transparent 
in the static-baryon sector (^71). However, for the ansatz made below, the matrix G is 
proportional to nnity: G = glNf ■ (Note that G transforms nnder the chiral transformation 
(pop as G I—*• UGU~^, so the mnltiples of nnity are hxed points of this gronp action.) If 
one decides to consider only those conhgnrations of Z and Z'^ for which G is scalar, the 
static-baryon action ( 0 ) simplihes to 

“Fbaryon [-^j “FYacuum [-^j S T COnst , 

SO the saddle-point expansion is rigoronsly jnstihed (for large Nc) if the integral is restricted 
to these conhgnrations. We will nse it to approximate the fnll integral. 


4 Vacuum saddle point configurations 

The saddle-point analysis for the action fnnctional S'vacuum[-Z^] has already been carried ont 
in [g, ^ 1 ^, so we are going to be brief here. In varying the action (|57|), the complex 
matrices Z and Z^ are to be considered as independent, which leads to two sets of eqnations. 
Variations of Z(n+|) ahect only the blocks M{n) and iV(n+|), with the linear response 
being 

5^vacuum = Tr {-M{n)-^ + Zt(n+f )iV(n+f )-') 6Z{n+§) . 

The resnlting saddle-point eqnation reads M{n)~^ = Z'^(n+|)iV(n+i)“^ or, by taking the 
inverse on both sides. 


M(n) = Z(n+i) + Zt(„+4)-‘ . 


(40) 
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Similarly, the variation 5Zl(rt-|) influences M{n) and N{n-^), and yields the saddle-point 
equation M{n)~^ = which is equivalent to 

M(n) = Zt(»-|) + Z(»-4)"‘- (41) 

As an immediate corollary, we have 

Z(„+|) + Zt(„+|)-1 = Zt(„-|) + (42) 

at every site n and for any pair /i, v. 


4.1 Homogeneous vacuum 

The simplest possibility for the held Z, is the scalar ansatz Z = = zl, with 2 ; a 

spacetime-independent real number. The vacuum saddle-point equations (|40|) and 
are solved by this ansatz if we put 


= 


1 + 


{amY 


am 


(43) 


If m > 0, the action S', 
powers of (am), we get 


vacuum 


\/2d + 1 V 2d + 1 2d + 1 

takes different values on these solutions. Expanding it in 


S„m[z±] = L^^^Nf 


In 


{2d + 2Y \ 

{2d + iY+yy ^ d + i 


V2dTl 

=F —^—:—am 


(44) 


which shows that for a positive quark mass, the configuration z = z_|_ minimizes the action. 

In the chiral limit (m = 0), a continuous set of solutions is obtained by applying the 
transformations (^^ to the homogeneous configuration Z = = ZvacI for 

Zvac = (2d+ 1)“^/^ . 


This vacuum configuration is invariant under the transformations of the diagonal subgroup 
Ui = U 2 E U(A^/) of the chiral symmetry group, but it maps to a new “vacuum” by taking 
Ui = U E U 2 = ^Nf By the Goldstone mechanism, the breaking of the continuous 

\]{Nf) symmetry leads to the existence of massless modes, namely the mesons, an effective 
Lagrangian for which was obtained by expanding S'vacuum near this vacuum in |]T^ 

If ±I 7 ^ 17 G U(A^/), the vacua obtained by translating the homogeneous one by U are 
staggered, in the sense that the value of Z depends on the parity of its position. However, 
on adopting the notations 17(n+|), lT(n'+§) for fields on the even and odd sublattices, the 
staggered vacua become homogeneous for each sublattice: 

!/(„+!) = {2d + 1)-V2[/ ^ lT(n'+f) = {2d + l)-B2t/t . ( 45 ) 


In [^, it was proved that, modulo the 'U{Nf) degeneracy, the configuration Z = ^(vac 
is the unique solution of the vacuum saddle-point equations in the chiral limit (except in 
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dimension d = 0, where the symmetry of the action is larger). The proof proceeds by a local 
argument, showing that for each site n the 2(d+ 1) saddle-point equations involving M{n) 
imply the equality of the matrices Z(n+|), Z^n-^) for all ju = 0,d; iteration of this 
result then trivially leads to the set of staggered configurations (^51). The proof strongly 
relies on being the Hermitian conjugate of Z, a constraint which is not mandatory. 
By relaxing it, we are now going to hnd a plethora of additional solutions of the vacuum 
saddle-point equations. 

4.2 Nonhomogeneous vacuum configurations 

By local considerations, as stated above, the only solutions of the vacuum saddle-point 
equations (and the Hermiticity constraint relating Z to Z^) in the chiral case are homoge¬ 
neous in both sublattice helds V and W. However, on a hnite lattice, say with the topology 
of a (d -|- l)-dimensional torus L'^ x T, there is also a global aspect to consider: one has 
to make a choice of boundary conditions for the various helds. The simplest choice are 
periodic boundary conditions in all directions, but one can also impose d-twisted boundary 
conditions, say along the hrst spatial direction 1: 

+ Li) = e'® H(n+§) , hh(n'+f + LI) = e-‘^ hh(n'+£) , (46) 

for all n and fi. An opposite twist for the helds V, W is natural in view of their opposite 
behavior under the chiral transformations (|30|). 

Now, accepting these twisted boundary conditions, let us investigate which conhgura- 
tion will minimize the action (|H^). A homogeneous conhguration suhers from a “phase 
jump” along a d-dimensional boundary, which is energetically very costly. A more rea¬ 
sonable ansatz for a minimum of the action is the following: the helds V, W smoothly 
rotate their phase, starting from V, hh ~ Zvac for ni = 1, to H Zyg^, W = Zvac at 
Ui = L, with a linear phase evolution in between. In this way, everywhere in spacetime 
the conhguration locally looks like one of the degenerate homogeneous vacua. 

4.2.1 Contour deformation 

The above ansatz for V, W is only qualitative. In order to actually obtain held conhgura- 
tions that satisfy both the twisted boundary conditions (^Bj) and the saddle-point equations 
(S,B. we need to relax the Hermiticity relation between the helds Z and Z^. 

By its construction via the color-havor transformation, the integrand is to be 

viewed primarily as a function of the real variables {{Zab + Zl)(.n+I), HZ^-Zl){n+l)}, 
the total number of which is D = 2Nj (d -I- 1) iTL*^). From |^, this function for 6 ^ 27rZ 
has no saddle points on M^, but it can be analytically continued into C^, where complex 
saddle points may exist. On such a saddle point there must exist at least one link (n,+|) 
where the matrix Z^ dihers from the Hermitian conjugate of Z. 

If a complex saddle point is not “too far” from the original contour of integration, 
it contributes to the vacuum-sector partition function, upon deforming the contour of 
integration so as to reach that point. 
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4.2.2 Vacuum saddle point equations for twisted fields 

We will demonstrate below the existence of complex saddle points for the vacuum sec¬ 
tor with any ^-twist. To make things simpler, we restrict ourselves to a 2-dimensional 
spacetime, with a twist in the spatial boundary conditions (we call the time index t, the 
spatial index x). The above qualitative ansatz for the fields V, W suggests the following 
symmetries: 

• All fields are scalar, i.e. at each point Z and are multiples of the identity matrix. 


• The fields V, W are time-independent. For each position x, there are 4 field variables 
associated with the time-like link which we denote by uo(a^), wo{x), Vq{x), Wq{x), 
and 4 field variables associated with the space direction, which we denote by Vi{x + 
1/2), wi{x + l/2), vl{x + ll2), wl{x + l/2). 


Thus, at each position a; = 0,... , L — 1 we have 8 independent complex variables. In the 
chiral limit (m = 0), the saddle-point equations (|40D pertaining to M{x,t) on an even site 
{x, t) read 


Vo{x) + Wq{x) + vi{x + 1/2) -f- - 1/2) = vo{x) + 1 /vq{x) 

= w*o{x)+ l/wo{x) 

= vi{x + l/2) + l/vl{x + l/2) ^ ’ 

= wl{x - 1/2) + l/'«;i(a; - 1/2) . 

The equations (^) pertaining to M{x, t-|-l) are obtained by interchanging u w, u* <-*• tc*. 
For a finite quark mass, 2am is to be added to the left-hand side. 


• The two first equations, together with their v w exchange analogs, allow us one 
more simplification. Indeed, they imply the identities Vo{x) = Wq{x), Vq{x) = Wo{x) 
(the alternative possibility, Uo(3^) = ^/wo{x) and Vq{x) = 1/wq{x), is incompatible 
with other relations that need to be satisfied). So there remain only 6 complex 
variables for each x. 


In the next section, we will provide approximate solutions for Eq. (^) together with their 
V w partners and assuming the above symmetries. 


4.2.3 Linearized problem 

To solve (at least approximately) the above equations, we will use the fact that we expect 
the fields to be locally close to one of the configurations (|^). We can then expand the 
saddle-point equations to first order in the perturbations from that configuration, and 
solve the linear problem. We start by expanding the fields around the real positive vacuum 


17 



V, W = l/^/^: 


Vq{x) 
Vi{x + 1/2) 
wi{x + 1/2) 



(1 + 5uo(a;)) , 

(1 + 5vi{x + 1 / 2 )) , 
(1 + 5wi{x + 1 / 2 )) , 


^;i(a: + l/2) 
wl{x + l/2) 



(1 + 5i;(;(x)) , 

(l + 5i;*(x + l/2)) 
(l + 5wt(a;+l/2)) 
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After inserting these expressions into (0) and expanding to linear order, we obtain a 
“transfer matrix representation” of these eqnations, i.e. a linear eqnation relating the vector 
of deviations of the spatial components of the helds {5ni, Sv^, 6wi, at position a;+1/2, 
to the same vector at position a: — 1/2. The strnctnre of the 4x4 transfer matrix allows 
to decompose it into two 2x2 matrices, npon considering at each point the vectors 


~ + Sw^J ’ ~ - 6w*J ' 

In terms of these two vectors, the linearized eqnations read 


(48) 


R(x + l/2) 
l(a: + l/2) 


(_( J) - 1/2) ='n R(i - 1/2) , (49) 

(l/2 “/f)Ki-l/2)='T,l(i-l/2). (50) 


Similarly, the deviations of the temporal components 6vo{x), Svq{x) are determined by the 
variations at a/ — 1/2: 


/ i5t'o(i)\ 

V4Ko(x)y 


/3/4 

f3/4 



■''■yy. -//.•)« 


1 / 2 ). 


(51) 


Thns, the transfer matrix allows to express the linear variation of all helds by the deviations 
at position 1/2. 

To make the az-dependence more explicit, we seek to diagonalize the transfer matrices 
Tr and Tj. The hrst transfer matrix T^ has the eigenvalnes 

— —(3 ± 2y/2), associated to the vectors R± ^ . (52) 

The second transfer matrix Tj cannot be diagonalized bnt only pnt in Jordan normal form. 
Indeed, it acts on the vectors 


def 




(53) 


18 


as 


Ti l+ = l+ + l-,T, L = L. 


Therefore, an initial deviation 


— C+r R+ + C-r R- ) Kl/2) — C+i 1+ + c_j l_ 


propagates throngh the transfer matrix as follows: 

R(x + 1/2) = c+ri- R+ + C_.(- e-^)" R_ , 
\{x + 1/2) = c+i 1+ + {c-i + X c+i) L . 


This linear evolntion is only valid as long as the deviations from 1 /\/3 are small compared 
to nnity. This cannot be the case nniformly for onr twisted ansatz, where the helds near 
X = L take valnes close to e=*"^® /v^. Still, the fact that l(x + 1/2) depends linearly on the 
position is enconraging: this is exactly the behavior we expect for the phases of the helds 
in the ansatz. 

The linearization of the saddle-point eqnations can actnally be performed near any 
of the degenerate family of vacna (^51). Linearizing the eqnations in the vicinity of a 
vacnnm v = e^^'^ /a/ 3, we obtain for the deviations the same transfer matrix as before. 
We can therefore constrnct local solntions near varions (p-vacna using (|5^, and glue them 
together to obtain a global, “rotating” solution. An equivalent procedure is to exponentiate 
the deviations. 


no (a:) 


1 

7 ! 


exp{5no(a;)} , 


(55) 


etc., and extend the equations for l(x + 1/2) to a larger domain of validity. This we 
do as follows. 

First of all, the R-part of the deviations grows exponentially, and is staggered with 
respect to x. Our ansatz excludes both features, so we simply set c+r = C-r = 0. Next, we 
note that the l-deviations depend on two coefficients, c+j and c_j. According to (^5l), their 
real parts describe the moduli of the helds, and the imaginary parts the phases. In our 
ansatz, we expect the moduli of the helds to be constant and close to l/\/3 (a linear growth 
would be incompatible with the boundary conditions). Therefore, we set Re c+j = 0. 

The other coefficient c_j causes a global shift of the helds, which can be interpreted 
as a “generalized chiral rotation”: the generalization consists in taking in Eqs. ( pUj ) for 
Ui and U 2 any invertible complex matrix, and replacing f/| and f/j by f/f ^ and U 2 ■ One 
easily checks that the ehective action is invariant under this GL{Nf) x GL{Nf) extension 
of the chiral symmetry group. The complex extension appears since we have relaxed the 
Hermiticity condition. The parameter c_i is then seen to parametrize the -manifold of 
scalar complex homogeneous vacua. 

The remaining coefficient, which we abbreviate to a *= Im c+j, is responsible for a 
nonhomogeneous solution. As we explained above, deviations l(x) ~ iax of order 0{1) can 
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be rescaled to deviations of order 0(a) by changing the reference vacunm e^^'^/A/S. As 
a resnlt, the helds bnilt from (^, with c±r = 0, Re c+i = 0 satisfy the saddle-point 
eqnations np to order (9(a^) nniformly in x. 

One can add terms of higher order in a to the exponent, so as to kill the higher-order 
terms in the expansion of (^71). By iterating the procednre, one obtains for the helds a 
series expansion in powers of a, snch that the saddle-point eqnations are satished to any 
order. We conjectnre that this expansion can be (re)snmmed, at least in a certain domain 
in a, thereby yielding an exact solntion of (^). The solntion np to order is 

' ' ^ „c_i-ia/2+5Q2/8 

Vi{x + 1/2) = ^ ^ 50 ) 

V3 

p—iax 

vUx + 1 /2) = e-c_i+ia-aV8 _ 

Vs 

The expressions for the helds ta, w* are obtained by replacing a —> —a, c_i —> —c_j. 

The coefficient a parametrizes the slope of the phase with respect to x, and it must be 
tuned according to the boundary conditions: 

aL = e + 2 'kQ^ , ( 57 ) 

where Qw is some integer. For a hnite twist 9, a can be chosen small only in the large- 
volume limit L ^ 1 , in which case a can take several values labeled by the integers 
< L. 


Vq\X) = 


73 


4.2.4 Topologically nontrivial configurations 


We now return to the original problem with periodic boundary conditions ( 6 ^ = 0). We have 
shown that there exist nontrivial solutions, for which the helds are position-dependent, 
with their phases rotating times when the position x goes from 0 to L. The integer 
can be called the winding number of the conhguration. We can associate a winding number 
to a (discrete) conhguration because the phases of the helds v, w are varying smoothly with 
position. More generally, when the lattice has the topology of a (1 + d)-dimensional torus, 
one can associate to any smooth scalar conhguration a set of winding numbers {Qw,ji}) 
each number specifying the number of times arg(n) rotates between the positions (n+|) 
and (n+L/i+l). 

Using Newton’s algorithm, we have searched for numerical solutions of the vacuum 
saddle-point equations (^7[), starting from trial conhgurations with Qyj = 1, Qyj = 2. We 
plot some results in Fig. These plots are very well described by our approximation 
including the (9 (q;^) corrections. 
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4.2.5 Nontrivial vacua for finite quark mass 

So far we have constructed nontrivial vacua only in the chiral limit, where a continuum of 
homogeneous vacuum configurations exists. What happens to these nontrivial vacua when 
the chiral symmetry is broken explicitly by switching on the quark mass? 

Recall that for m ^ 0 there remain only two homogeneous scalar vacua (^) out of 
the former U(l) continuum, with one of them {z^ac = ^+) being an absolute minimum of 
the action. One can again study the linearized saddle-point equations near this solution. 
Unlike before, the results will apply only locally, since we cannot use the trick of rescaling 
the homogeneous reference vacuum any more. 

We have performed numerical searches for (1 + 1)-dimensional topologically nontrivial 
vacua, assuming the same symmetries as before (we took scalar, time-independent, smooth¬ 
ly varying fields). The saddle-point equations are modified by the addition of the quark 
mass term 2am to the left-hand side of Eqs. (^). The outcome of these calculations (see 
below) can be understood in large part by analytical reasoning, as follows. 

To linearize the saddle-point equations around the point z+, we set 

vo{x) = z+ etc. (58) 


As in the chiral case, the 4x4 transfer matrix splits into two 2x2 matrices that apply to 
the vectors R, I defined in Eq. (^51). These matrices can be written for an arbitrary value 


of am, using the exact expression (^) for z+(am) (we only consider the case d = 1). Using 
the same notations as above, they are 


Tr{am) 

Tj(am) 


1 

1-4 

1 


f - z +^ - 2 - 4 24 + 24 \ 

V -24-24 -4 + 24 + 34^) ’ 
fz +^ -2 + 4 24 - 24 \ 

V-24 + 24 4 + 24-34^) ’ 


(59) 


and the deviations on links pointing in the time direction propagate as 


/ 5no(4^ 


2(1 - 4) Vl 



;f)K-i/2). (60) 


One easily checks that both transfer matrices have the property DetTj,/j(am) = 1. 

We expand both matrices and their spectra in powers of am, since we are interested in 
the case of a small quark mass. In the massless limit the matrix Tr(0) is hyperbolic, with 
negative real eigenvalues that are well separated from each other (one expanding, the other 
contracting). Thus, a perturbation of order 0{am) is still diagonalizable, with eigenvalues 
and eigenspaces shifted by that same order 0{am). We will keep calling the eigenvalues 
— e^^, with the expansion 


A(am) = ln(3 + 2\/2) + am/VG + 0{{am)'^) . 


(61) 


On the other hand, Tj(0) was nondiagonalizable with eigenvalue +1, so a perturbation 
can change its qualitative features. For any positive am, Tj(am) becomes diagonalizable. 
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with real positive eigenvalues associated to eigenvectors (g±^pam))- To express the 

deviations of the held vo(x), we use the coefficients 

drf 1 + 4 

1 + 4 

For small am, these data have the following expansions: 

7(am) = ^(amy^^ + ^(arn)T2 + 0((amy/y , 

Ki(am) = - ^(arn)T2 + 0((amy/y , (62) 

K±o(am) = + ^am ± ^ + 0((amy) . 

For a hnite mass am, one expects the linear approximation to be valid only for small 
deviations. However, the error introduced in the saddle-point equations by the linear 
approximation is at most of order 0{am\6z\), so it remains small if the mass is small. 

The numerical solution we obtained for a mass am = 0.01 and winding number = 1 
(see Figs. suggest that the vectors R are negligible in a large domain of x around the 
point X = 0 where helds are close to This indicates that the coefficients c±r vanish, 
like in the chiral case. The helds will therefore depend on two complex parameters e±: 

6vi{x + 1/2) = e+ eT'^ +e_ e”"^^ , 

5vl{x + 1/2) = -e+ e"i+^" -e_ e-'i-T'" , (63) 

6vo{x + 1) = e+ e''+°+^^ +e_ . 

As opposed to the chiral case, both coefficients e± are complex, so that both the phases 
and the moduli of the helds vary with x. This ansatz hts the numerical solution even when 
the deviations from become of order 0{1), which is quite surprising. However, it is 
unable to reproduce the zone where the helds cross the negative real axis (near x = L/2). 
For a quark mass am = 0.01 the values of the various exponents are 

7 = 0.1527 , = -0.3045 , k+q = -0.0673 , k_o = 0.0845 . (64) 

These values are used in the hts to the numerical solution shown in Fig. 


5 Static baryon saddle point equations 

The ehective action Fbaryon for the static-baryon sector, Eq. (|37|) , contains a “string” term 
in addition to the “sea” term S'vacuum- While the sea term depends on every one of the 
matrices Z(n+^), the string term involves only those matrices Z and that are situated 
in the near vicinity of the string. More precisely, what enters into the baryon world line 
propagator, G, are the matrices iV((i+i)6) and M{to). Of these, the former depend only 
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on Z and along the string, whereas the latter also involve the matrices Z(to+/i/ 2 ) and 
Zl(to-/i/ 2 ). The Z held on the remaining links (away from the string) appears only in 
‘S'vacuum, SO the Variation with respect to these matrices yields the same equations (^) and 
(PI) as in the vacuum sector. 

For simplicity let us consider the two particular cases N^. = 2 and Nc = 3, using the 
expressions (^) for the effective action. The most general variation yields 


Nr = 2: = 


baryon — ^*^vacuum 


7V, = 3: = 


baryon — ^*^vacuum 


Tr(G5G) + Tr(G)Tr(5G) 

Tr(G2) + (TrG)2 ’ 

(TrG)2Tr5G + 2TrGTr(G5G) + Tr(G2)Tr(5G) + 2Tr(G25G) 
(TrG)3 + 3TrG Tr(G2) + 2TrG3 


We then work out how the various traces of powers of G respond to variations of each 
matrix Z and entering in the dehnition of G. For instance, variations of Z{td+fi/ 2 ) with 
/i 7 ^ 0 affect only the matrix M{td), whereas varying Z((t+i/ 2 )o) affects both M{td) and 
A^((t+i/ 2 ) 6 ). These computations are simplihed by the use of cyclicity properties: given 
any decomposition G = G 1 G 2 , we may replace G in the static-baryon action by the matrix 
G = G 2 G 1 , as G always appears under a trace. 

We provide detailed calculations for the variation with respect to C = Z{{t+i/ 2 )d). The 
modihed factors of G in this case are and A^((t+i/ 2 )o), and the modihed matrix G 

reads 


G + 5G = ■ ■ ■ + 5C{- + C^N{{t+i/ 2 )d)-^)}N{{t+i/ 2 ) 6 ) ■■■ . 

It is now natural to conjugate G + 5G into 

r(G + (5G)T-i = G + SC{ - N{it+i/ 2 )d)-^)G , 

where G iV(p+i/ 2 ) 6 )M((t+i) 6 )“^ ■ ■ ■ iV((t-i/ 2 ) 6 )M(t 6 )“^. The saddle-point equation that 
follows from varying Z{{t+i/ 2 ) 6 ) then takes the succinct form 

6 Z{{t+i/ 2 ) 6 ) : 0 = ( - M(f6)“^ + Z'^((t+i/2)6)A^((t+i/2)6)“^) ■ (ijv^ - F)v^(G)) , 

with the case-by-case dehnition of the matrix-valued function F)v^(G) being 

" ■ ^(^ 2 ) + (TrG )2 ’ 

AT . Fin - ^(TrG)2 + 2G'2TrG + GTiiG^) + 2G^ 

(TrG)3 + 3TrGTr(G2)+ 2Tr(G'3) ‘ 

The saddle-point equation obtained by varying Zl((t-i/ 2 ) 6 ) is similar. It is best expressed 
in terms of the matrix G' = M(t 6 )“^iV(p+i/ 2 ) 6 ) ■ ■ ■ M(p-i) 6 )“^iV((t-i/ 2 ) 6 ): 

0 = (Iat^ - F)v^(G')) • ( - + A^((i-i/2)6)“^Z(p-i/2)6)) , 
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with the same dehnitions for F/v^ as above. The term multiplying the unit matrix 
stems from ^S'vacuum- Note that this term factors out in both of the above variations. 

To get an idea of the matrix ^Nf — F]s[^{G), we compute it in the vacuum conhguration 
Z = Zva.J-Nf In this case we have G = G = G' oc Itv^. We then notice that F^^^alNf) = 
Nj^ ItVj. for any number of colors and any a ^ 0. Therefore, on the vacuum conhguration, 
we get Iatj. — Fn^{G) = {1 — which is invertible as soon as Nf > 1. More generally, 

this equation holds as long as G is a multiple of the unit matrix, which is a property of the 
inhomogeneous scalar ansatz we will make in the next section. 

Clearly, as long as the matrix — F/v^(G) remains non-singular, the saddle-point 
equations due to varying Z((i+i/ 2 )o) and Zt((t+i/ 2 )o) are identical to those in the vacuum 
sector, Eqs. (^) and (H). As was said earlier, this is also the case for the equations due 
to varying all matrices Z and not involved in the matrix G, i.e. those away from the 
string. The only difference to the vacuum equations comes from the matrices on the links 
adjacent to the string, namely Z{td+fi/ 2 ), Z^{t 6 -ji/ 2 ) for the directions /x = 1,... ,d. These 
matrices are contained only in some M~^ factor of G, and their variations give the following 
saddle-point equations: 


6 Z{t 6 +fi/ 2 ) : —M{t 6 ) ^ + Z\td+ji/ 2 )N{t 6 +fi/ 2 ) — M{tb) ^Fjv^(G), (65) 

5Zl(t6-A/2) : + A^-^(i6-A/2) Zitb-i^li) = - Fjv,(G') , (66) 


where the matrices G, G' are the same as before. These equations represent the only 
obstruction that prevents the vacuum conhguration Z = from being also a 

saddle point of the static-baryon sector. 


5.1 Configurations for the static baryon in 1 + 1 dimensions 

In this section we present approximate solutions of the saddle-point equations in the static 
baryon sector for the simplest nontrivial case, which is the two-dimensional Euclidean 
square lattice {d = 1). We use the same notations and assume the same symmetries as in 
Section 4.2.2| , so at each position x there are 8 independent complex scalar variables. 

The baryonic string is placed on the Euclidean time axis at position x = t). The 
equations to solve are the vacuum saddle-point equations (^) off the string (x 7 ^ 0 ), and 
the mo diked equations 


Uo(0) -f wS(0) + Ui(l/2) + wt(-l/2) + 2am = uo(0) + l/n*(0) 

= m;( 0 ) + l/mo( 0 ) 

= (i-Af7‘)A.(i/2) + i/«:(i/2)} 

= (1 - Af7‘){»;(-i/2) + iM(-i/2)} 

on the string, together with the equations obtained by exchanging v w, v* ^ w*. 
As in the vacuum sector, these equations imply the identideations vo{x) = Wq{x) and 
Wo{x) = Vq{x) for all x. 
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5.1.1 Physical requirements 


Recall from Section |4.1| that demanding to be the Hermitian conjugate of Z (and 
assuming the vacuum saddle-point equations) leads to a homogeneous configuration, where 
the fields are constant on each sublattice. Such a homogeneous configuration cannot satisfy 
the last two of the equations (|67D . To get a solution, we must relax the Hermiticity condition 
(cf. Section |4.2.1|) , and consider the helds Z and Z'^ as independent variables. 

We want the baryon to be a localized object, in the sense that a baryonic saddle-point 
conhguration should differ from a vacuum configuration only in some neighborhood of the 
baryon world line. The baryon can then be interpreted as a spatially localized excitation of 
this vacuum. A priori, baryon excitations may exist on top of each of the vacua described 
in Section |. 

In Eqs. (^), the number of flavors Nf enters just as a parameter, so one can extend 
the equations to any real value oi Nf. In the limit Nf = oo, we recover the vacuum 
saddle-point equations. We can therefore obtain a solution of the baryon saddle-point 
equations by starting from a given vacuum conhguration (at Nf = cx)), and deforming 
the conhguration by continuous variation oi Nf down to its physical value (say Nf = 2). 
Any baryon conhguration obtained in this way carries the same topological charge as the 
vacuum it is associated to. 


5.1.2 Topologically trivial sector = 0, chiral limit 

In the sector with zero winding number, we numerically found a unique solution (see Fig. 
asymptotic to the homogeneous vacuum ^^vac, he. satisfying the asymptotic condition 

* * Ixl—>c» ^ , /— 

W’ = 1/V3 . 

All the helds of this conhguration are real and time-independent. Various components 
coincide pairwise or in quadruples: 

Vq = Wq = Vl= Wl= Zq , Vi = Wi = Zi, Vl=wl = zl . 

One easily checks that the saddle-point equations are invariant under the following trans¬ 
formation: 


Zi{x + 1/2) zl{-x-112) , zo{x) ^ zo{-x) , ( 68 ) 

which represents a rehection at the baryon world line. The solution found numerically 
is invariant under this transformation, and we believe the same to be true for the exact 
solution (or else we would get a second solution by rehection). We can therefore restrict 
our study to the domain 0 < x. 

Our numerics show an exponential convergence of all the helds towards z^s.c as we 
depart from the string (see Fig. and the signs of the deviations alternate with x. This 
phenomenon can be explained by the linearized saddle-point equations studied in Section 
[4.2.3| . The linear theory indeed applies if the helds are close to Zvac, which is the case for 
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large enough x. Eqs. (^, ^), together with the physical condition that the deviations 
decay as a: —>• cx), require c+r = 0. We thus get the ansatz 


zo{x) — ® ’ 


Zl{x + 1/2) — + C-r (- 


^—X\X\ 


(69) 


.-:(x+i/2) = ^(i+c_. 


e^ (— e , 


with e^ = 3 + 2\/2 as before. According to this linear approximation, the helds oscillate 
around the asymptotic value Zvac, and the amplitude of the oscillations is controlled by a 
unique coefficient, which we denote by Ci = 


The results (|6^) £t the numerical conhguration not only far from the string (where this 
is expected), but even down to the baryon string, where the fields deviate significantly 
from z^sjc- More precisely, the ansatz fits zi, for all x, whereas zq departs from it only at 

X = 0. The value 2:o(0) *= ^ (1 + (Tq j together with the parameter Ci can be computed 


using the (nonlinear) saddle-point equations on the string (|67|) and the reflection symmetry 
(|68|). We obtain two equations: 

Cq + 2C'iC'o + 2Ci + 4:Cq = 0 , 

3e^ + 4e^ CqCi + 4Co + Ci{3 + 7e^) + 4 = 0 . 

The equations have four pairs of solutions, two real ones and two complex ones, conjugate 
to one another. The physical solution (which deforms to Ci = Cq = 0 as we vary Nf from 
2 to cxo) is Cl = —0.0971 and Cq = 0.0504, giving 


V2e^Ci = -0.8002 . 


(70) 


The relative smallness of these constants (except for the last one, which governs the am¬ 
plitude (5zo(0)) may explain why the ansatz works well down to the string. 

The above saddle-point configuration is indeed situated outside of the original contour 
of integration: it is a “complex saddle point” (although all fields have real values). A 
contour deformation has to be performed for the variables i{zi — zl){x), which move away 
from the real axis as |x| decreases. 

This ansatz is tailored to the limit L —>• cx), but owing to the fast decrease towards Zy^c, 
it is already quite good for short lattices. In Fig. | (bottom), we show a logarithmic plot 
of the deviations of the helds from the homogeneous vacuum, for a lattice of total length 
L = 20, as well as the values predicted by the above ansatz. 

This conhguration might be called a non-topological soliton, cf. 


211. Its characteristic 


length (in units of the lattice spacing a) is = 0.5673, and its mass will be computed in 
the next section. 

This real scalar conhguration is just one point on a U(77/)-manifold of solutions, ob¬ 
tained by the action (^OD of the chiral symmetry group. For a generic point on this manifold. 
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the configuration is staggered in time. These solutions are saddle-point configurations for 
Nf = 2 , and do not depend on the number of colors Nc (which does not appear in the 
saddle-point equations). However, the value of the action for these configurations does 
depend on see the next section. 


5.1.3 Topological baryon, chiral limit 


We also obtained topologically nontrivial configurations, characterized by a nonvanishing 
winding number Q^^,. In Fig. ^ we plot a solution in the baryon sector with = 1. All 
fields are scalar, and have the symmetries described in Section [4.2.2. 


Away from the string (x 1 ), the moduli of the fields are close to Zvac ^ind the phase 
varies linearly. In this region, we can apply the linear theory described in Section [4.2.3| , 
in particular the ansatz (|^). Now only the coefficient c+r has to vanish, to prevent the 
deviations from exploding as x —> oo. As in the vacuum case, the coefficient c_j can 
take any value, yielding only a chiral shift of the fields. We find that c+j = ia is purely 
imaginary, as in the vacuum. The fields are well fitted by 


vi{x + 1 / 2 ) = + c_r (- 

V3 




etc. 


( 71 ) 


for positive x. Near the string, the moduli of the fields behave in a similar way as in the 
non-topological sector, whereas their phases make a small jump at x = 0. Both this 
jump and the value of c_j, can be computed from the full saddle-point equations near the 
string (see below). The value of the parameter a depends on the height of this jump: a 
will not be exactly equal to its value in the vacuum, which is 211 jL for this topological 
sector, but rather to (27r — (3)/L. As a consequence, the convergence of the fields towards 
the vacuum configuration away from the string will not be exponential, but only linear 
(the fields coincide at the “antipode” of the baryon, x = T/ 2 ). 

From Fig. |^, we can assume that the above ansatz is still a good approximation for 
ni(l/2), and for v\{l/2) up to a phase jump of {3/2. Given this assumption and setting 
Nf = 2 , the saddle-point equation on the baryon worldline, no( 0 ) + l/nQ( 0 ) = l/ 2 {ni(l/ 2 ) + 
l/nj‘(l/ 2 )}, yields two real equations, one of which reads 


sin (3a/2 +/3/2) 3(1 + Gi) 

sin(a/ 2 ) 2(1 + Gie'^) 


For small angles a and /3, this equation gives a linear relation between them: using the 
value for Ci obtained in the non-topological sector, we get [3 ~ 3.24 a. The value of the 
slope a then is a = 2iiQu,/ (3.24 + L). 


5.1.4 Topological baryon, m 7 ^ 0 

As in the vacuum sector, our results for a finite quark mass are mostly numerical (see 
Figs. We obtained solutions of the saddle-point equations with various winding 

numbers, which are close to the corresponding vacuum configurations except in the vicinity 
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of the baryon string. Near the antipode of the string (x = L/2), the held approaches the 
homogeneous value z+, and the helds can be htted by the linear theory developed in Section 

OH - 


6 Mass of the static baryon {Nf = 2 ) 

The mass Mbaryon of the baryon is dehned by comparing the static baryon partition func¬ 
tion to the vacuum one. Since the saddle-point conhgurations are classihed according 
their winding number the comparison is performed within a given topological class, 
i.e. between a baryon conhguration and the corresponding vacuum sector. 

In the limit of large lattices, the ratio of partition functions is expected to behave as 


^baryon,Q,. {L X T) ^ j. 

'2^vacuum,(5tj; (LxT) 


OO 


(72) 


As explained in Section p.2| , we estimate both partition functions through their respective 
lowest-order saddle-point approximations: 




baryon 




exp 


( -^c{*S^baryon['^baryon] *S^vacuum['^vacuum]}) • 


(73) 


As before, we will only treat the (1 -|- l)-dimensional case. 


6.1 Mass of the static non—topological baryon 

We hrst study the baryonic excitation on top of a homogeneous vacuum, which has van¬ 
ishing winding number. The value of the action for the homogeneous vacuum was given in 
Eq. It is proportional to the volume of the lattice, and is called a “sea” term in the 

literature. 

The action S'vacuum is part of the full static baryon action ([37|), which leads to a “sea” 
contribution to the baryon mass: 


M def “Svacuum ['^baryon] ‘S'vacuum ['^vacl 

sea = Nc -^^^ 


We now assume the limit m = 0 (Section ^.1.2 ). By the time-independence of both 
conhgurations, the sea contribution reads 


L/2 


M,ea = Ar^iV, I ^ In 

:X = —L/2 


N,{x + l/2)No{x) 
M{x) 


— In 


^—^baryon 


iVi(x + l/2)iVo(x) 
M(x) 


Z=Zv 


Since the deviations of the helds from Zvac are small and decrease exponentially — see 
Eq. (|BUD — it is reasonable to keep only the linear order, and extend the sums to L = 
OO. Quite remarkably, this linear approximation gives a vanishing sea term (for m = 
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0). Alternatively, we can compute the sea term from the Nf = 2 conhguration obtained 
numerically. In this way we obtain 

Msea ~ 0.02324 X Nc (in units of a“^) . 


This answer is small (5%) compared to the second term we compute below, so the linear 
approximation (giving Mgga = 0 ) is rather good in this respect. 

The remaining contribution to the baryon mass comes from the sum over traces of G, 
which involve only the helds on or adjacent to the string. In the QCD context this is 
generally referred to as the “valence quark contribution” |21|. For the time-independent 


conhguration described in Section |5.1.2|, the valence term for Nf = 2 becomes 


M, 


valence 


= In 


( ao /Nc + Nf — 1 


I cq 


N. 


+ Nc In 


2zo(0) + 2zi(l/2) \ 

l + ;^o( 0)2 J 


^ha 


(74) 


The term proportional to Nc evaluates to 


In + In 


1 + C'i/2 + Cq/2 
1 + Co/2 + Cq/A 


In 73 + Ci/2 = 0.5493 - 0.0485 


We notice that the second term due to the deviations l^^baryon — ^vac| is small compared to 
the hrst term ln\/3, obtained by inserting the vacuum conhguration ; 2 vac into (|74D . 

The result of this linear approximation is very close to the exact (numerical) value, 
0.5004. On including the combinatorial and normalization terms (see Eq. (^)), we hnally 
get, for the non-topological static baryon in the chiral limit: 


Afbaryon ~ 0.5236 Nc + ln(l + Nc/2) + ln(l + Nc) ■ 
This yields for example 


(75) 


Nc = 2 : Mbaryon = 2.839 a ^ , 

Nc = 3: Mbaryon = 3.873 a~^ , 

where we have reinstated the mass scale given by the lattice constant. 


6.2 Masses of topological baryons 

To compute the baryon masses in the topologically nontrivial sectors, we hrst need to eval¬ 
uate S'vacuum ou the Corresponding vacua with winding number Qn,. This is straightforward 
in the chiral limit, where we have the accurate approximation (^6|) at our disposal. The 
result up to second order in the small parameter a = 27iQyj/L is 


T (*Fva,r,niim I 


•^vacuum I ^ vac, a 


‘S'vacuum ['2^vac]) -7/"^ 


3q(^ 


AL 


-Nf 


29 












The vacuum energies are obtained from this by multiplication with the number of colors, 
Nc- They agree well with the values computed numerically for L = 120 (and Nc = 2): 

= 0.25, = 0.98. 

The baryon mass in each topological sector is dehned relative to the corresponding 
vacuum energy Eq\ For a nonvanishing quark mass, both the vacuum energy and the 
baryon energy are computed numerically. In the following table, we summarize our results 
for a lattice L = 120, with Nc = Nf = 2, for various quark masses: 


Qw 

m = 0 

am = 0.002 

am = 0.01 

0 

2.839 

2.840 

2.852 

1 

2.840 

3.097 

3.692 


In Fig. im we plot the baryon masses in the topological sectors = 0 and Qw = 1 as a 
function of the quark mass. 

7 Zigzag baryon 

In the standard formulation of the theory on the Euclidean 2-dimensional square lattice, 
the temporal and spatial directions are given by the lattice generators, t = ( 1 , 0), a; = (0, 1 ). 
However, it is also possible to use different spacetime axes. For instance, on the same square 
lattice, we define the zigzag spacetime axes as f = (1/2, —1/2), x = (1/2,1/2) (see Fig. 0 ). 
The unit spacetime separation now has the length a/\/2. For convenience, we choose the 
spacetime origin on the middle of a link (say, a link in the direction t + x), so that the 
coordinates of lattice sites will be half-integers, while links will be indexed by integers. 

The division of the square lattice into two sublattices is now expressed only in terms 
of the spatial coordinate: the links at even positions x = 2n carry the helds V, 1/1, while 
the links at odd x = 2n + 1 carry the helds IF, IFh 

The vacuum effective action and its corresponding saddle-point equations are still given 
by the formulas after a suitable change of labels for links and sites. 

7.1 Time—independent vacua 

Once again, we make a scalar time-independent ansatz for the helds: 

V a: even, V t : V{x,t) = v{x)I , V\x,t) = v*{x)l , (76) 

V a; odd, V t : W{x,f) = w{x)I , W\x,t) = w*{x)l . (77) 

In particular, this implies the equality of helds situated on links (a;, f) and (a;, t + 1). There 
are 2 scalar variables at each position. 
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(78) 


With these symmetries, the vacuum saddle-point equations read: 

X even : 2v{x) + 2w*{x — 1) + 2am = v{x) + l/v*{x) 

= w*{x — 1 ) + l/w{x — 1 ) , 

X odd : 2w{x) + 2v*{x — 1) + 2am = w{x) + l/w*{x) 

= v*{x — 1) + l/v{x — 1) . 

In the chiral limit, the homogeneous vacuum conhgurations are still given hj v = w* = 
6 “^ 2 ;vac and v* = w = z^nc- For helds close to this vacuum {v{x) = e“^ exp{5u(a;)} 

etc.) the linearized saddle-point equations yield the following transfer matrix equation: 


6z{x) 

dz*(x' 


f-3/2 -1/2 

I 1/2 - 1/2 


= T. 


Sz{x — 1) 


(79) 


6z{x — 1) 

6z*(x — l)J ^^^\5z*(x — 1)^ 

The symbol 6z stands for either Sv or 6w, depending of the parity of x. In contrast with 
Section [4.2.3| , we now have just one transfer matrix, which relates deviations of v,v* to 
deviations of w,w* and vice versa. This transfer matrix Tzig is related to the matrix Tj 
described in Section |4.2.31 Indeed, it acts on the vectors I+, L as follows: 

'1^ I- 1 To-o^l— I- . 


zig^ + 


The deviations ^u(O), 5u*(0) are parametrized by two complex parameters c± as 

Gt'(o))=<=+'++<=-1- 

The deviations will then depend on position as follows {x is even): 


6v{x) 

6v*(x] 


c_ + (x + l)c+ 
—c_ — {x — l)c+ 


6w{x + 1 ) 
6w*{x + 1 ) 


-c_ - (x + 2 )c+ 

C_ + XC+ 


( 80 ) 


The rest of the discussion is identical to the one following Eq. (|55|) . The coefficient c_ 
plays the role of a global shift, or “generalized chiral rotation”. If Re c+ ^ 0, the absolute 
values of the helds vary linearly with x, which is incompatible with their periodicity. On 
the other hand, taking c+ = ia will linearly rotate the phases of the helds, keeping them 
close to some vacuum conhguration, as in Eqs. (^). Taking for a a multiple of 27r/L, we 
obtain a topologically nontrivial conhguration. 

The case of broken chiral symmetry (m ^ 0) can be treated along the same lines as 
in Section |4.2.5| ; the above linear evolution in position is then replaced by an exponential 
one, at least for helds in the vicinity of the value 2 ;+. The transfer matrix takes the form 
/2 - 1/2 


„-2 


(- 


1/2 
1 

exp{±K.zig) 


- 34/2 


, and has the eigenvalues — exp(± 72 jg) associated to the eigenvectors 


), with the expansions 


'^zig 


^zig 


V 2 

31/4 

2 V 2 

■ 3 ^ 


amY^'^ + 0{{amY^Y , 
amY^’^ + O^iamY^Y . 


(81) 

(82) 
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Up to this order, the exponent '^zig differs from the corresponding exponent 7 of Eq. (| 6 ^ 
by a factor of a/2. This factor actually compensates for the ratio of unit lengths between 
the two frameworks (a versus a/\/2). For the quark mass am = 0.005 used in compiling 
the figures, we have ^zig = 0.0762 and e'^^'s = 0.85905. 


7.2 Zigzag baryon 

In the new labeling conventions, the worldline of a static baryon situated at position 
a; = 0 forms a zigzag curve (see Fig. O). Assuming that the fields are scalar and time- 
independent, the G matrix appearing in the baryonic part of the action is 

G = {(2v(0) + 2w*{-l) + 2am)-/l + n(0)n*(0)) 

(2w(l) + 2t;*(0) + 2am)“/l + n(0>*(0))}^/2 I 

The resulting saddle-point equations on the string read 

2n(0) + 2 m;*(-1) + 2am = 'u(O) + l/n*(0) 

= (i-]v/‘){w-(-i) + iM-i)}, 

2m(l) + 2!)*(0) + 2om = !i*(0) + 1 /k(0) 

= (l-]V7‘){w(l) + l/u,-(l)}. 

The saddle-point equations are invariant under the transformation 

w{—x) < —>■ w*{x) , v{—x) < —^ v*{x) , (84) 


which is also a symmetry of our numerical solutions. 

In the chiral limit, the linear dependence in (|80|) makes it impossible for the absolute 
values of the fields to approach z^ac at infinity, unless the deviations are purely imaginary; 
this latter possibility (/(a;)| = z^ac) is incompatible with the saddle-point equations on the 
baryon string (PB]). However, for a finite lattice, infinity is the “antipodal point” x^o = T/2. 
In numerical searches (see Fig. [^, top), we found a solution which comes close to z^ac near 
the antipode (but does not converge exponentially to it). For x even, 0 < x < L/2, the 
fields are well described by 


v{x) 
v*{x) 
w{x + 1 ) 
w*{x + 1 ) 


— ^vac 
'^vac 
'^vac 
'^vac 


^c+(x+l-L/ 2 ) 

q-c^{x-1-L/2) 

^-c+{x-\- 2 -L/ 2 ) 

qC+{x-L/2) 


(85) 


where the value of the real coefficient c+ is small (for a lattice of length L = 80, we found 
c+ ~ 0.04). The fields in the region —L/2 < x < 0 are obtained by applying the symmetry 
(|^). At all points x / 0, the fields are close to a “generalized homogeneous vacuum”. 
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The coefficient c+ is fixed by using the above ansatz up to x = — 1, and then enforcing the 
equations (^): in this way one obtains (to lowest order in c+) the transcendental equation 

c+ e^+^ = 1 , (86) 

with approximate solution c+ ~ L~^ InL. The field on the baryon takes the value 

v(0) = v-{0)^^^. ( 87 ) 

This configuration is very “distorted” compared to the original contour of integration. 
Indeed, the ratios v/v* and w/w* are of order L in the vicinity of the string. In the case 
L = 80, the above equations yield c+ = 0.04018 and n(0) = 0.1736, in excellent agreement 
with our numerical data. 

7.2.1 Zigzag baryon, broken chiral symmetry 

If chiral symmetry is broken by a nonvanishing quark mass, the deviations from 2 ;+ evolve 
exponentially with x' = x + L/2 away from the antipodal point x' = 0. More precisely, we 
should have, in some domain |x'| <C T/2, 

/ fc(x') \ ^ / e+(—+ e_(— 

y5x*(x')/ y—e_,_ e'^^®s(— Q-i^zig(^— Q-'izig^x' 

From the mirror symmetry (^dj), the coefficients e± are related by 

e_ = e+ . (89) 

We numerically computed a solution with L = 80, am = 0.005 (see Fig. |^, bottom), for 
which this ansatz works well up to the string. Using this fact, it is possible to estimate the 
value of e+ (the last remaining parameter), as in the chiral case. The crudest approximation 
yields the equation 

-e+ K^igC 0 ^*^"+ = 2 , with C = exp{'y^ig{L/2 - 1)} . 

For our configuration, the solution of this equation is e+ = 0.0614. This configuration has 
a mass M = 2.885 a~^ with respect to the homogeneous vacuum. 

There also exist topologically nontrivial vacuum and baryon saddle-point configurations 
with the symmetry (^4|). For example. Fig. |T^ shows the solution in the sector = 1 with 
quark mass am = 0.005 on a lattice of length L = 160. With respect to the corresponding 
vacuum configuration, this configuration has a mass M = 3.18 a~^. 

8 Concluding remarks 

The “color-flavor transformation” introduced in |T^, ^ replaces an integral over the gauge 
group U(W) by an integral over the “flavor” degrees of freedom. In the present paper we 
extended this transformation to the gauge group SU(iVc)- 
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The color-flavor transformation can be interpreted as a kind of duality, linking two dif¬ 
ferent formulations of the theory. We believe that this duality transformation may be useful 
for treating realistic non-perturbative QCD. Here we have applied it to a simple model of 
two-dimensional lattice fermions. The non-Abelian theory we have treated is of course too 
far from realistic four-dimensional QCD for our results to be of direct phenomenological 
relevance. 

The main approximation we made was to assume the strong-coupling limit for the 
lattice gauge fields. In this approximation gluons do not propagate, and the connection to 
asymptotic freedom at short distances is lost. An unphysical consequence is the absence 
of the U(l) chiral anomaly. Thus, the chiral symmetry group of our low-energy effective 
action is not SU(A^/) but the larger group U(A^/). 

Among the gauge groups SU(Ai"c) the case Nc = 2 is special, as the vector and covector 
representations of SU(2) happen to be equivalent. Since these representations correspond 
to quarks and antiquarks respectively, there is no physical distinction between baryons 
and mesons in that case. This symmetry between baryons and mesons is obscured in the 
present treatment which, by the use of a saddle-point approximation valid only for Nc'^ 1, 
is geared to the large-A'^c limit. It can, however, be made manifest by identifying SU(2) 
with the compact symplectic gauge group Sp(2) and using the color-flavor transformation 
for the latter 
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To extend the formalism to lattice QCD in four dimensions, we need to take into account 
the spin degrees of freedom and put the chiral fermions properly on the lattice. We hope to 
address these issues in a separate publication. Here we only note that a first step towards 
a more realistic color-flavor transformed theory of the strong interaction was described in 
P^ , Q , where we discuss the effect of spontaneous chiral symmetry breaking and estimate 
the numerical values of the chiral condensate, the pion decay constant and the mass of the 
pion. 
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A Action of the color and flavor groups 

Transformations U G GL(A'"c) of color space act on the one-fermion operators by 

= fUU-'f . , 

TuflJu" = fUuf , = !UU-^f . 

Transformations e GL(2Nf) of flavor space can be decomposed in the way shown 

in Eq. (|^, and the action of the various factors may be described separately. An element 
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c 


1 z 
0 1 


exp Q j acts by 


y rk rri—l rk Ik ry 

-^CJ-\-b'^(^ J-\-b J—a^ba : 

TATy 

rn Ik rp — \ 7k 

-^CJ+b-^C — J+b 7 

Tj^Ty 


an element diag(y4, D) 


e GL(2]V/) by 


fl, + fiA . 

fk 

J-b 5 


rn rk rp—\ _ 'T' ^—1 _ £k T\ 

dia.g{A,D)J +b-‘- dia.g(A,D) ~ J+a\^ Jba ’ diag(A,D) Jdiag(A,D) “ J-a^ab ^ 

rn £k rp—\ _ Ik A rji Jk rp—\ _ 7k { r) — ^\ 

ddiag(yl,D)J+6-i diag(yl,D) ~ J+a^ab , diag(A,D) J-fed diag(A,D) ~ Jba ’ 

1 O' 


and an element = 


Z 1 


e GL{2Nf) by 


rn_ rk rji — l _ rk 

~ J+b 5 

T i_ Ik rji—l _ 7k I £k ^ 

(^J +6"^ ^ J -\-b J—a'^o.b •) 


k rji — 1 

c 

k rji — 1 

c 


T^f-bT^ 

T^f-bT: 


rk 

J-b 5 

Ik rk y 

J—b J • 


All these formnlas are particular cases of the fermionic Fock-space representation of the 
Lie group GL(2iV/) expounded in Chapter 9 of [|T^]. 


B Normalization constants 


We are going to calculate the normalization constants = J^dg KEglTglEQ)]"^ intro¬ 
duced in Eq. o - for the values Q = 1, 0, — 1. To that end, we employ the decomposition 
of the group G = U(2iVj) given by Eqs. ([T3|) and (p^. This yields 


\{Bi\Tg\E,)f 


\{Vl + ZZ^),aUai\^^^ 

Det(l + 


for Q = 1, and similar expressions for the other two cases. The first step now is to do the 
integral over U G U(A^/), which for Q = ±1 is effectively an integral over a {2Nf — 1)~ 
dimensional sphere. Carrying it out by the method of Section p.5| , we get the preliminary 
expressions 


, f dZdZi 

° ' Jc"!-', Det(l + ZZt)2«,+». • 

(N, - l)!iV,! f [{1 + Z'Z)nf‘ dZdZ' 

^’(N^ + Nf-iy. Det(l + ZZt) 2 y,+y. 


(90) 

(91) 


where Cnj is defined by 


^ M. — 


dZdZ^ 


'N. 


f Det(l + ZZf)^^f 


(92) 
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For later convenience, we have made a change of integration variables Z <->■ ZMn the 
numerator of the integral in (|9T|). 

In the second step we perform the integration over the Nf x Nf matrix Z using a 
recursion procedure similar to that in |^. From here on we use the simplihed notation 
n = Nf. The recursion consists in slicing the matrix Z into vertical vectors, step by step. 
We now detail the hrst step of the recursion. We decompose Z as Z = (Zn^n-i, ^i), where 
Zi is a (column) n-vector, and is a n x (n — 1) matrix. We then have the expressions 

(Z<Z)n = . 

Using the (positive dehnite) n x n matrix Fi which is dehned as the square root of 

F^ = 1 + Zn,n-lZl n-l ) 

we make a change of variables, from zi to Wi = From 1 + ZZ'^ = Fi(l + tcir(;|)Fi, 

we get the relation 

Det(l + ZZ^) = (1 + wjwi) Det(l + Zn,n-iZi „_i) . 


The change of variables from Z to {Zn,n-i,wi] has the Jacobian Det(l + Zn,n-iZ\.^_i). 
Each of the integrals (|^) , (|^) , and (^) can now be written as the product of a Zn,n-i- 
integral times a tci-integral. 

The former can in turn be expressed as the product of a „_ 2 -integral times a W 2 - 
integral (with W 2 a n-vector), which can be decomposed in turn, and so on, until we reach, 
at the n-th step, a i-integral, i.e. an integral over the hrst column of the original matrix 
Z. We call this column vector Wn for reasons of homogeneity. 

The successive Jacobians multiply to give the following integration measure: 

dZdZ"^ = dw\dwi (1 + W 2 W 2 )dW 2 dW 2 ■ ■ ■ (1 + wl^WnY~^dWndwl^ . 


The integrands in (|^p^ also have simple expressions in the new variables, due to the 
identities (Z'^Z)ii = wl^Wn and 

Det(l + ZZ^) = (1 + wlwi){l + W 2 W 2 ) ■ ■ ■ (1 + wlwn) . 


The tCjUntegrals to be performed are all of the type {N > n) 

f dw^dw „{N — n)\ 


Icn (1 + 


= TT 


m 


The resulting expressions for the normalization constants are 

1 {2Nf + N,-l)\---{Nf + N,)\ 

1 Nf{2Nf + Nc-l)\---{Nf + Nc + l)\ 

it "? (iV/-l)!'"0! ’ 


(93) 

(94) 

(95) 
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where we have reinstated n = Nf. The quantity entering into the baryon mass is the ratio 


Nf 

ao Nf + Nc 


(96) 


C static baryon 

In this appendix we prove the formula ([TtI) for the action functional of the static baryon 
sector. We need to integrate polynomials in the quark helds along the world line of the 
baryon (the baryon “string”), weighted by the same Gaussian as in the vacuum sector. 


We start out using the short-hand notation fO = 0 + fO of Section p.l.2| for sites and 
links on the string. The part of the integrand containing the quark helds situated on the 
string, namely V’(tO) for f = 0,... , T — 1, then reads 


Xi(V’(0), iV(i6)^/>(10)) Xi(V’(10), Ar(f6)V’(20)) ■ • • 

■ • • Xi(V’(T-b 6 ), ^(T- 5 ) 6 )'?/'(T 0 )) X exp — ^ V;*(f0)M(f0)'?/i*(t0) . 


T-l 


t=o 


Isolating the terms with fermions at the site n = tO, we are faced with the integral 

f d'4){n)d^p{n) Xi{'4’{n— Q),N{n-b/ 2 )ip{n)^ (^^(ri)^ iV(n+ 6 / 2 )-^(n + 6 )) 


= tti 


(]v!?]v/-l )!) j n<'’/4(n)#i(") e 


-bcG)A^cc' 
(*)/ 


X ^ sgnasgnr - b)Nah{n-b/ 2 )i)l^'’\n)\Y'^i'iN)^a'h'{n+b/ 2 )'il)lr\n + ^) 

o-,Te6iv^ i j 

2 

sgn(ar) JJ “ 6 )Ar„.fe^(n-o/ 2 ) 


= 


(7V^-1)! 




CTjTSSjVc 


di}\n)di)\n) |iV„/fe/(n+6/2)i/>7^(n + 0) 

where the hrst equality sign uses the expression (|2T|) for the function xi- Note that the 
integral between curly brackets involves only fermions of color i. The fermionic version of 
Wick’s theorem yields for it the value Mr^,{n) Det M(n), so after combining the permuta- 
tions a and r, the above expression becomes 

(N, - 1)P 


a 


\N, + Nf-l)\^ 

Nc + Nf -1 


= ai Det M(n) 


NJ Det MY sgn p ^jJl. (n - (n-b^n+b)^p^!;"^ (n + 0) 

pe&Nc * 

-1 

Xi (V’(n — 0), G{n-d->n+b)7p{n + 0)) , 


N- 


f 
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with the “propagator” G(n-o^n+o) N{n-b/2)M{n)~^N{n+o/2). 

Repeating the procedure, we successively integrate over the quark fields along the string, 
by which process the matrices N and M get organized into a single propagator. In the 
final integration step, we need to take into account the periodic boundary conditions for 
the quark fields: ipiTO) = '0(0). The final integral over 0(0) then reads 

J d0(O)d0(O)xi(0(O),G(o^T6)0(O)) . 


We now use the following expression for the function xv 


Xi(0,0) = ai 




Nc 


E n 


(Nc + Nf - 1 )! ^ 

' / / a&eMc 




which is easily obtained from Eq. (|^) by interchanging the product over colors with the 
sum over flavors. Wick’s theorem then yields for the '0(O)-integral the result 


0.1 


(iVc + iV^-l)! 


DetM(0)'^^(-l) 


\N, 


{ai,bi}a£&Nc * 




bAo- 




The last matrix product may also be expressed in terms of the propagator G defined in 
Eq. (ID, G = G(o-To)M(0)-h 

What’s the interpretation of the sign factor (—1)'^'"? To answer that question, recall 
that we evaluated the Grassmann field integral using time-periodic boundary conditions 
(instead of the conventional t ime-ant ip er iodic ones). In a d-dimensional quantum mechan¬ 
ical frame work with Hamiltonian H and inverse temperature (3, this would mean that we 
are computing not the usual partition function but rather the supertiace Tr (—l)^^e“^'^ 
with Np the total fermion number. The overall sign factor (—1)-^= originates from that 
very fermion number, and is simply telling us that the baryon is a fermion (boson) if W is 
odd (resp. even). 

Let us take a closer look at the contributions from the sum over permutations a G 
Each permutation a can be uniquely decomposed into a product of independent cycles. 
Denoting by q(ct) the number of cycles of length I in this decomposition, the contribution 
from a to the partition function can be written as 


Nc 


Nc 




Ci(o-) 


{ai} 2=1 


1=1 


The permutation group &Nc be partitioned into disjoint classes with respect to con¬ 
jugation (cr, a' are said to be conjugate to each other iff there exists a permutation r such 
that a' = r“^crr). Two permutations a and a' are in the same conjugacy class iff they 
have the same cycle structure, i.e. V/ : q((j) = ci{a'). This allows to rewrite the sum over 
cr as a sum over the conjugacy classes a E &Ncj taking into account the cardinality of each 
class, AA(d), given in Eq. (|^). We then obtain the result (|37D . 
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Figure 3: Vacuum configuration computed numerically, with winding number = 1 in 
the chiral limit. Fields which are numerically indistinguishable (e.g. |ni| and |tCi|) are 
represented by the same symbol. There is a perfect £t with formulas (^), including the 
correction. The difference between |no| and |ng| comes from Re c_i ^ 0. 
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Figure 4: Numerical vacuum configuration with winding number = 1 and chiral sym¬ 
metry explicitly broken by a hnite quark mass. The helds are normalized with respect to 
the corresponding value z+. The phases arg(no) and — arg(nQ) are indistinguishable, so we 
plot them together (idem for arg(ni) and — arg(t(;i), resp. arg(t(;j|‘) and — arg(nj‘)). 






Vacuum am=0.01 Qw=1 (L=120) 
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Figure 5: Same vacuum configuration as in the previous figure, plotted on a logarithmic 
scale. We fit the fields in the range 0 < |a:| < 30 (where they are close to z^) with the 
ansatz (|63|), using the theoretical values for 7, ki, from Eqs. (^). The best fit is 
obtained with e+ = (—2.80 + i7.62) x 10“® and e_ = (—1.82 — i5.48) x 10“^. 








Z/Zvac 


Baryon m=0 Qw=0 (L=120) 



Baryon m=0 Qw =0 (L=20) 


deviations (log-plot) 



Figure 6: Numerical baryon configuration in the chiral limit. All fields are real. Top: the 
fields converge exponentially fast to z^ac- Bottom: we compare the numerical data (circles, 
triangles, squares) with the theory of Section b.l.2| (3 lines, cross for 2:o(0)). 
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Figure 7: Numerical baryon in the chiral limit with winding number = 1 (we only 
plot the vicinity of the baryon worldline). The absolute values are very similar to the case 
= 0, but now the phases vary linearly away from the worldline. The slope a and the 
phase jump (3 at the baryon are in good agreement with the theory of Section ^.1.3 
hnd 27r/a = 123.42 and argno(l) — argno(—1) = 0.046 5 . 24 q !. 
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Baryon am=0.01 Qw=1 (L=120) 
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Figure 8 : Numerical baryon with winding number = 1 and broken chiral symmetry. 
Away from the worldline, the fields converge exponentially fast to The logarithms 
of the phases are linear in the range 3 < a: < 55, with slopes 0.1523 < 7 < 0.1533 in 
agreement with the theoretical value from Eqs. (| 6 ^) . 

















Figure 9: Numerical baryon with winding number = 2 
The symmetry with respect to a; = 0 is only approximat 
accuracy. 
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Figure 12: Numerical zigzag baryon for the case with chiral symmetry (top) and without 
(bottom). We plot the (real) helds on a logarithmic scale. Top: a linear £t yields the slope 
c+ = 0.04 and the held on the baryon n(0) = 0.1735, in excellent agreement with formulas 
(1^,^). Bottom: we use the exponential ansatz of Section 
over the domain |a;'| < 30. The values e+ = 0.062, e_ = —0.053 are in good agreement 
with the analytical theory. 


7.2.1, with coefficients e± fitted 
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Figure 13: Numerical zigzag baryon with winding number Q.^ = 1 and broken chiral sym¬ 
metry. We plot the real (top) and imaginary (bottom) parts of \og{z/z^) on a logarithmic 
scale, together with a £t by the ansatz (|88D in the region \x'\ < 60. The coefficients take 
the values e+ (3.48 -|- il.07) x 10“^, e_ —(3.00 -|- iO.92) x 10“^, and satisfy quite well 
the relation (p9|). 











